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Abstract 

We present the two-loop virtual QED corrections to e + e~ — > and Bhabha scattering 

in dimensional regularization. The results are expressed in terms of polylogarithms. The form of 
the infrared divergences agrees with previous expectations. These results are a crucial ingredient 
in the complete next-to-next-to-leading order QED corrections to these processes. A future 
application will be to reduce theoretical uncertainties associated with luminosity measurements 
at e+e~ colliders. The calculation also tests methods that may be applied to analogous QCD 
processes. 
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1 Introduction 



Bhabha scattering is an important process for extracting physics from experiments at electron- 
positron colliders primarily because it provides an effective means for determining luminosity. These 
measurements depend on having precise theoretical predictions for the Bhabha scattering cross 
sections. As yet, the complete next-to-next-to-leading order (NNLO) QED corrections needed for 
reducing theoretical uncertainties have not been computed. In this paper we present the complete 
two-loop matrix elements that would enter into such a computation. This calculation also provides 
a means for validating techniques that can be applied to physically important but more intricate 
QCD calculations. It also provides an additional explicit verification of a general formula due 
to Catani Q for the structure of two-loop infrared divergences, and allows us to determine the 
process-dependent terms for the processes at hand. 

In Bhabha scattering there are two distinct kinematic regions: small angle Bhabha scattering 
(SABS), and large angle (LABS). In the LEP/SLC energy range, SABS is used to measure the 
machine luminosity via a dedicated small angle luminosity detector. SABS has a large cross section 
— about four times larger than Z decay in the 1° — 3° window — making it particularly effective as 
a luminosity monitor. At the same time, SABS is calculable theoretically with high accuracy from 
known physics (mainly QED), apart from hadronic vacuum polarization corrections that rely upon 
the experimental data for e + e~ annihilation into hadrons at low energy p| |||. Therefore, SABS is 
an important ingredient in measuring any absolute cross section. For instance, the measurement of 
the hadronic cross section at the Z peak, <rg, which enters several precision observables, is especially 
dependent on an accurate theoretical understanding of Bhabha scattering. 

At LEP/SLC, large angle Bhabha scattering interferes with e + e~ —* Z — > e + e~ and so it is 
needed to disentangle important parameters such as the electroweak mixing angle. It is also useful 
for measuring the luminosity at flavor factories such as BABAR, BELLE, DA<£NE, VEPP-2M, and 
BEPC/BES [B]. A peculiarity of future electron linear colliders is that the luminosity spectrum 
is not monochromatic because of the beam-beam effect. Because of this, measuring the total 
small angle cross section of Bhabha scattering alone is not sufficient, and therefore the angular 
distribution of LABS was proposed for disentangling the luminosity spectrum M. 

Due to the experimental importance of this process, significant effort has been devoted to 
developing Monte Carlo event generators — see for instance ref. Q for an overview. In order to 
match the impressive experimental precision, a complete inclusion of NNLO QED quantum effects 
has become necessary. On the theoretical side, however, the calculation of two-loop four-point 
amplitudes has been a roadblock to further progress. 

In this article we present the two-loop virtual QED corrections to the differential cross section for 
Bhabha scattering, i.e., the two- loop amplitude interfered with the tree amplitude and summed over 
all spins. We neglect the small electron mass in comparison to all other kinematic invariants, and use 
dimensional regularization to handle the ensuing infrared divergences. Besides these contributions, 
a number of other virtual and real emission contributions (discussed in the conclusions) still need 
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to be obtained before a full Monte Carlo program for the Bhabha scattering cross section can be 
constructed. 

The two-loop QED four-fermion amplitudes are also a useful testing ground for two-loop QCD 
calculations containing more than one kinematic invariant, which are required for higher-order 
jet cross sections and other aspects of collider physics. For processes that depend on a single 
momentum invariant, a number of important quantities have been calculated up to four loops, 
such as the total cross section for e + e~ annihilation into hadrons and the QCD /3-function @. In 
contrast, the only complete two-loop four-point scattering amplitudes presently known for generic 
kinematics in massless gauge theory are the N = 4 super- Yang- Mills amplitudes ||, [|, and gg — > gg 
in a single helicity configuration in pure gauge theory |in]. The two- loop amplitudes required for 
NNLO computations of jet production in hadron colliders, or for NNLO three-jet rates and other 
event shape variables at e + e~ colliders, remain uncalculated. We note in passing that partial results 
for the leading-color part of two-loop contributions to quark-quark scattering have very recently 
appeared [11]. 

Two important technical breakthroughs are the calculations of the dimensionally regularized 
scalar double box integrals with planar [^] and non-planar pj| topologies and all external legs 
massless, and the development of reduction algorithms for the same types of integrals with loop 
momenta in the numerator (tensor integrals) [14, 15, 16, 17, 11]. Related integrals, which also arise 
in the reduction procedure, have been computed in refs. [IS, 19]. Taken together, these results are 
sufficient to compute all loop integrals required for 2 — ► 2 massless scattering amplitudes at two 
loops, thus removing a major obstacle to several types of NNLO calculations. In this paper we use 
these techniques to evaluate the integrals encountered in the Bhabha calculation. An even more 
recent result concerning two- loop planar double box integrals with one massive external leg [20] 
holds promise for the NNLO computation of three-jet rates at e + e~ colliders. 

There has also been significant progress in developing general formalisms for other aspects of 
NNLO computations involving massless particles. The motivation has typically been infrared-safe 
observables in QCD, but many of the developments can be applied to the Bhabha process as well. 
The developments include an understanding of the intricate structure of the infrared singularities 
that arise when more than one particle is unresolved (i.e., is soft or collinear with another parti- 
cle) pi] , |22| , |23| ]. Improved approximations to the NNLO correction to splitting functions have been 
constructed recently as well [24|. 

Infrared divergences are a significant complication in all the QCD and QED computations men- 
tioned above. In any suitably "infrared-safe" observable all final-state divergences will cancel [25]. 
However, divergences occur in individual amplitudes for fixed particle number, and it is very useful 
to have a general description of such divergences. Catani has presented a general formula for the in- 
frared divergence appearing in any two-loop QCD amplitude Q. By appropriately adjusting group 
theory factors, it is straightforward to convert Catani's QCD formula to a QED formula, allowing 
us to directly verify it. Moreover, we extract the exact form of a process-dependent term in the 
formula, for the case of QED scattering of four charged fermions. Previously, the only process for 
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which this term had been extracted was the quark form factor which enters Drell-Yan produc- 
tion [26]. (It should also now be possible to extract it for gg — > Higgs using the recent two-loop 
computation |2^].) Interestingly, a simple generalization of the quark form factor term (converted 
to QED) correctly predicts the process-dependent term for the e + e~ — > fi + and Bhabha am- 
plitudes. We also use Catani's formula to conveniently organize the infrared divergences and to 
absorb some of the finite terms. 

The previously computed non-abelian gauge theory amplitudes [^, ||, 10] were obtained via 
cutting methods. The low multiplicity and relative simplicity of the e + e~ — > ^ fj," and Bhabha 
scattering Feynman diagrams makes it relatively easy to directly compute the diagrams, as we do 
here. We include here only the pure QED diagrams, neglecting for example the contributions of 
Z exchange, and hadronic vacuum polarization effects. The former are negligible at this order in 
SABS and in LABS at flavor factories. The hadronic contributions are important, but much of 
their effect is straightforward to include by introducing a running coupling. 

We perform the calculation in dimensional regularization [28] with d = 4 — 2e and set the 
small electron mass to zero, since it is the only form in which the required two-loop momentum 
integrals are known. Moreover, it provides a powerful method for simultaneously dealing with both 
the infrared and ultraviolet divergences encountered in gauge theories. Traditionally, dimensional 
regularization is not used for QED, in part because the infrared divergences are relatively tame 
compared to non-abelian gauge theories, so photon and electron masses are sufficient for cutting off 
the theory. Another important reason for using dimensional regularization is to validate techniques 
that can also be applied to the more complicated case of QCD. In QCD, dimensional regularization 
is the universally utilized method for dealing with divergences. 

In the high-energy Bhabha process, even with an "infrared-safe" (calorimetric) final-state def- 
inition, the electron mass will still appear in large logarithms of the form L = ln(Q 2 /ml) due to 
initial-state radiation. However, in the dimensionally regulated amplitudes these singularities (like 
all others) appear as poles in e. It may therefore be most convenient to handle the initial-state 
singularities using an electron structure function method [29| implemented in the MS collinear 
factorization scheme. 

In the next section we briefly describe our method for computing the two- loop amplitudes. 
Then we describe Catani's formula for the divergence structure of the amplitudes, followed by a 
presentation of the finite (C(e )) terms for both e + e~ — ► fJ. + n~ and Bhabha scattering. In the 
final section we give our conclusions, including some discussion of the remaining ingredients still 
required for construction of a numerical program for Bhabha scattering at this order. 



2 The Two-Loop Amplitudes 

The 16 independent Feynman diagram topologies describing the two-loop QED corrections to 
e + e _ — > and Bhabha scattering are enumerated in fig. |l|. In this figure we have suppressed 

the fermion arrows. After including the fermion arrows and distinct labels for the external legs, 
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there are a total of 47 Feynman diagrams; however, many of these diagrams generate identical re- 
sults. Of the 47 diagrams, 35 contain no fermion loop, 11 contain one fermion loop, and 1 contains 
two fermion loops. The Bhabha amplitude may be obtained from the e + e~ — > /i + jjT amplitude 
by adding to it the same set of diagrams, but with an exchange of one pair of external legs. The 
e~ ' fjT — > e~/i - and e~fi + — » e~ fi + amplitudes may, of course, be obtained by crossing. 









Figure 1: The independent diagrammatic topologies for two- loop four- fermion scattering in QED. 



We have evaluated these diagrams interfered with the tree amplitudes and summed over spins 
in the conventional dimensional regularization (CDR) scheme. This interference gives directly the 
twodoop virtual correction to the 2—5-2 differential cross section. The rules for implementing CDR 
are straightforward because all particle are treated uniformly in all parts of the calculation. In this 
scheme, all momenta and all Lorentz indices are taken to be D = 4 — 2e dimensional vectors. (The 
7-matrices remain as 4 x 4 matrices; i.e., Tr[l] = 4.) 

After performing all 7-matrix algebra present in the twodoop Feynman diagrams, we use the 
conservation of momenta flowing on the internal lines to express the tensor structure of the diagrams 
in terms of inverse scalar propagators and a small number of additional scalar invariants containing 
loop momenta. The inverse scalar propagators cancel propagators in the denominator to generate 
simpler "boundary" integrals. To handle the integrals containing scalar invariants, we introduce 
Feynman parameters and interpret the resulting integrals in terms of scalar integrals with multiple 
propagators, which are then reduced to a set of master integrals with the help of equations in 
refs. gj, m, Hi. 



Proceeding in this way, we obtain an expression for the amplitude in terms of master integrals 



(of the type listed in ref. [15], plus a few more for the planar double box topology) multiplied by 
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coefficient functions. This expression is in principle valid to an arbitrary order in e, assuming that 
the master integrals could be evaluated to such an order. However, it is a bit too lengthy to present 
here, and for NNLO computations only the series expansion in e through O(e ) is required. To carry 
out this expansion, we use expansions of the master integrals presented in refs. [12, 13|, 14. [l5|. 
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As noted in ref. [11], there is a slight problem with the original choice of basis [|14|] for the two master 
planar double box integrals. In that basis, the coefficients for generic tensor integrals contain 1/e 
poles, necessitating an O(e) evaluation of the master integrals. Several solutions to this problem 
have been presented fll7| , |Tl|[ . We have used a slightly different solution, which is simply to use 
the original pair of master integrals defined in ref. [14|, except evaluated in d = 6 — 2e instead of 
d = 4 — 2e. In d = 6 — 2e the integrals have neither ultraviolet nor infrared divergences, making 
them simpler to evaluate through O(e ) than the d = 4 — 2e integrals. 

Many of the master integral expansions quoted in refs. 0, ||, @, ||, ||, || are in terms of 
Nielsen functions |3C], 



hi 



n—1 



t ln p (l - xt) , 



(2.1) 



(n-l)!p! Jo t 

with n + p < 4. We have found it useful to express the results instead in terms of a minimal set of 
poly logarithms [31|, 



Li n (x) 
Li 2 {x) 

with n = 2, 3, 4, using relations such as 
Si3(x) = 



x f x dt 
Jo t 



(2.2) 



Li 4 (l -x) + ln(l - x) Li 3 (l - x) + - ln 2 (l - x) \Li 2 {x) 



( 2 ) + ^ln 3 (l -x)lnx + C4, 



S22OE) = Li 4 (x) - Li 4 (l - x) + Li 4 



1 



ln(l 



Li3(a:) - 
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for < x < 1. 



+ -i- ln 4 (l - x) - - ln 3 (l - x) Inx + -C2 ln 2 (l - x) + ( A , 



Here 



c 2 



n=l 



7T 



1.202057. 



Ci 



7T 

90 



(2.3) 
(2.4) 



The analytic properties of the non-planar double box integrals are somewhat intricate [|13|| , since 
they are not real in any of the three kinematic channels for the 2 — > 2 process, 



s-channel : 
t-channel : 
n-channel : 



s > 
t > 
n > 



t, u < , 
s, n < . 
s, i < 0, 



(2.5) 
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where s = [k\ + k^) 2 ^ t = (k\ — k^) 2 , and u = (ki — k%) 2 . Therefore we shall present explicit 
formulae for the finite terms in the amplitude in both the s- and u-channels; those in the t-channel 
will be related by symmetries. 



2.1 General Structure of Divergences 



Dimensionally regulated two-loop amplitudes for four massless fermions contain poles in e = (4 — 
d)/2 up to 1/e 4 . The structure of most of these singularities has already been exposed by Catani Q, 
who described the infrared behavior of general two-loop QCD processes. We shall therefore adopt 
his notation in presenting our results. 

We work with ultraviolet renormalized amplitudes, and employ the MS running coupling for 
QED, a(p?). Of course this scheme can always be converted to another one, for example a(Q 2 ) 
defined via the photon propagator at momentum transfer Q, by a finite renormalization. The 
relation between the bare coupling a u and a(^ 2 ) through two-loop order can be expressed as PJ 



a u iXq S e = a(fj?) [i 



2c 



2\ A) . 2 1 2\ (Pi A 



2c 



+ 0(«V)) 



(2.6) 



where S e = exp[e(ln47r + ^(1))] and 7 = 
coefficients of the QED beta function are 



>(1) 



0.5772 ... is Euler's constant. The first two 



N f 
3it 



N f 



(2.7) 



where Nf is the number of light (massless) charge 1 fermions. 
The renormalized four-fermion amplitude is expanded as 

-M 4 (a(/x 2 ),/i 2 ; {p}) = 4vra(/i 2 



2tt 

^)VV;M)+0(<*V)) 
The infrared divergences of a renormalized two-loop amplitude in QCD or QED are Q, 

|M 2 V;M))r.s = i (1 W;M) |M 1) (M 2 ;M))Ra 

+ /^( e ,/i 2 ;M) |M V;{p})>rs + lM 2)fin (^ 2 ;M))aa , 



(2i 



(2.9) 



where \M^n\n 2 ] {p})) R .a is a color space vector representing the renormalized L loop amplitude. 
The subscript rs stands for the choice of renormalization scheme, and \x is the renormalization 
scale. These color space vectors give the amplitudes via, 



M n (l ai ,...,n a ") = (a 1 ,...,a n \Mn{pi,...,p n )) 



(2.10) 



where the are color indices. The divergences of M n are encoded in the color operators (e, fi 2 ; {p}) 
and I^(e, // 2 ; {p}). In the QED case, the color space language is clearly unnecessary; M.^ and 
are just numbers. 
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In QCD, the operator I^\e, [i 2 ; {p}) is given by 

J (1 W 2 ;M) 



1 e -^W » " 

EE t <- t j 



2r(l-e) 



e 2 + T? e 



M 2 e -iAy7rx 6 



2p t • pj 



(2.11) 



i=l j'^i 

where A*j = +1 if i and j are both incoming or outgoing partons and Aj,- = otherwise. The color 
charge Tj = {T"} is a vector with respect to the generator label a, and an SU(N C ) matrix with 
respect to the color indices of the outgoing parton i. The values required for QCD are 



T2 rri2 s~1 

g — 1 g — 0_P ■ 



T 2 g = C A = N c , 
lg = l ~lc A - 2 -T R N f 



(2.12) 



For QED we let Ca — * 0, Cp — > 1, Tr — > 1 and Tj • Tj — » e^ej = ±1, where the ej are the 
electric charges, to obtain 



/W(e jM 2 ;M) 



-6^(1) 



rfi 



2 3 

^7 + - 



2 \ e 



/i 



2 \ e 



+ 



/i 



2 \ e 



-It 



for the four-fermion amplitude 

e + (ki) e~(k 2 ) -» /u + (/c 4 ) n~(k 3 ) . 

(Note that the charges of incoming states should be reversed in computing Tj • Tj. 
The operator 1^.1 is given in QCD by Q 

/L 2 W; M) = -^ {1) (^ 2 ; W) (/ (1 W 2 ;M) + ^) 



(2.13) 



(2.14) 



+ 



e +^(i)r(l - 2e) ^2^/?o 
HS(e^ 2 ;{p}) J 



+ K J( 1 )(2e, j u 2 ;{p}) 



(2.15) 



where the coefficient is: 



if 



'67 _^ 
18 6 



C A - ^T R N f 



(2.16) 



For the QED process ( |2.14 ), we insert T 1 ) from eq. ( 2.13| ), take (5q and f}\ from eq. ( |2.7| ), and let 
K -» -107V//9. 

(2) 

The function -ET^a is process-dependent but has only single poles: 

ffS(e,/x 2 ;W) = o(iA). 



(2.17) 



Ref. |lj] does not give an expression for JTrs for a general amplitude, but only for the case of a 
qq pair, i.e. a single charged fermion pair. The result, which is extracted from the two-loop QCD 
computation of the electromagnetic form factor of the quark 2£], is 



H 



(2) 

qq,CUR 



2e r 



4er(l-e)l 2 Pl - P2 



1 56 

- 7(1) + 3C F if + 5C 2 vtA)Cf - — vr^ C F 



16 



7(3 CfCa 



(2.18) 
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where 



(1 7 88 \ / 4 32 \ 

-y - yC 2 + 2A( 3 )C F C A +(- + —( 2 )C F T R N f . (2.19) 

Performing the usual conversion to QED yields a result applicable to the electromagnetic form 
factor of the electron, 



H 



(2) 



e+e - , CDR 



1 e- 6 ^ ( /i 2 e" iAl27r 
2pi ■ p 2 



2e r 



4eT(l -e) 



+ 6C 2 - 12C 3 + 



25 
27 



+ C 2 UVf 



(2.20) 



Using our two-loop computation, and an all-orders-in-e computation of the one-loop amplitude 
for e + e~ — > (see sect. |2.2[ ), we have verified that the singular behavior of the e + e~ — > [x + 

amplitude in CDR agrees precisely with that predicted by eq. ( |2.9D in all three kinematic channels.^ 
In addition, we have extracted the function H^ e _ , CDR controlling the 1/e poles in eq. ( |2.9| ). 
We obtain 

r(2) 



e+e _ /i+/i _ ,cdr. 



1 e^ 1 ) 
4eT(l-e) 



x 2 



S 2, +(C 



+ 6C2 - 12C 3 + 



2 \ 2e 



25 

27 



(2.21) 



+ C2 N f 



This result agrees with a "naive" generalization from the form factor case, in which one sums 



eq. ( 2.20 ) over the six pairs of charged legs in the four-point amplitude, weighted by the sign of 
the charge product e^ej. (Note that the factors of (// 2 /(— Sij)) 2e are purely conventional here, since 
their deviation from unity only contributes at the level of finite parts, (D(e ). However, the overall 
normalization is predicted correctly by the sum over the six pairs.) 



2.2 e + e — > fi + fi at One Loop to All Orders in e 

In order to verify the structure of the infrared singularities, and to extract the finite remainder 
of the two-loop amplitude presented below, we computed the one-loop e + e~ — > fi + pT amplitude 
(interfered with the tree amplitude) to all orders in e. The result is 



£MM 0) M-£M 0) M 0)t + 



spins 



spins 



£=1 



(2.22) 



where the first term is the MS counterterm, expressed in terms of the tree-level interference 



EM 0) M 0)t = 8 



spins 



t 2 +u 2 



(2.23) 



1 Strictly speaking, we have computed the interference of the two-loop amplitude with the tree-amplitude, summed 
over intermediate fermion spins, so in our verification eq. (2.9) should be similarly understood to be interfered with 
the tree amplitude. 
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and 



A« = 4£(1 - 2eH (2 - 3e)u 2 - 6etu + 3(2 - e)t 2 Box^(s,t) 
s z L J 



1 - 2e s 2 



(4 - 12e + 7e 2 )i 2 - 6e(l - 2e)tu + (4 - lOe + 5e 2 )« 2 Tri(t) 



2\„,2 



1 



(1 - 2e)(3 - 2e) s 



2e(l - e) i((l - e)i - eu) N f - e(3 - 2e)(2 - e + 2e 2 )t« 



+ (1 - e)(3 - 2e)(2 - (1 - £)e + 2e 2 )i 2 



Tri(s) . 



The symmetry operation S acts as 



(2.24) 



S : t <-> u, £ 



(2.25) 



After carrying out the operation of S 1 , one should then set £ = 1. (Basically, ^ allows us to separate 
diagrams based on whether they have an even or odd number of photons attached to the muon 
line. Because photons have C = — 1, this criterion governs the t <-> u symmetry properties.) 

In eq. (|2724|) , Bax (6) (s, t) and Tri(s) are one-loop box and triangle integrals, the former evaluated 
in an expansion around d = 6 — 2e. For the divergence formula ( [2.9D , we need their series expansions 
in e through 0(e 2 ). In the it-channel, where the functions are manifestly real, their expansions are 
given by 



Box( 6 )(s,t) 



-l-e 



2e 2 (lu(-v) + WLi 3 (-v) - \v 2 U 2 (-v) - - \v z W + \v 2 W 2 

t2 -2 



6 



Tri(s) 



(-«)• 



1-e 



\V 2 -^VW-2C^J + (*<->t) 



12 3 16 



where 



s 1 t/ i f s 
w = — , u; = — , l/=m 

u u \ u 



W = In 



(2.26) 



(2.27) 



The expansions in the s- and t-channels can be found using analytic continuation formulae such 
as p| 



ln(l — x + ie) = ln(x — 1) + in , 

LI2 (x + is) = —L12 ( —] In 2 x H h ivr In x . 

\x J 2 3 



1 



Lia(x + ie) = Li3 I — I In 3 x H lnx + i— In 2 x , 

\x J 6 3 2 

/ 1 \ 1 7T 2 7T 

Li 4 (x + ie) = -Li 4 1 - ] - — In 4 x + — In 2 x + 2(4 + i- In 3 x , 
\ x / 24 6 6 

x > 1 , 



(2.28) 
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where ie is an imaginary infinitesimal added to s,t or u before continuing. 

We have verified that through 0(e°) our result for the one-loop amplitude agrees with a pre- 
vious calculation p^], up to terms which can be identified as being due to the conversion between 
dimensional regularization and a photon mass regularization. 

2.3 Modifications for Bhabha Scattering 

In comparison with the process e + e~ — » [i + fj>~ described above, the Bhabha scattering process 

e + (fci)e-(fc2) e + (k 4 )e-(k 3 ), (2.29) 

has additional exchange diagrams. In general, the interference required for Bhabha scattering is 
given by 



m{ Li) m 



(Li) A/ ,(X 2 )t 



spins 



Bhabha 



(la) K ~a(Li)\ 



spins 

+ u 



spins 



EM M, M ,T +EM ^ 

spins spins 



(Li) k ~ a (L 2 )] 



(2.30) 



where the symmetry U acts as 



U : s <-» t, 



(2.31) 

is the L-loop amplitude for e + e~ — > fi + fi~, and is the same L-loop amplitude but with 

legs 1 and 3 interchanged (taking into account the Fermi statistics minus sign). 

2.4 Bhabha Scattering at One Loop to All Orders in e 

In the CDR scheme, the tree- level exchange contribution required for Bhabha scattering in eq. ( 2.30) ) 
is 

Y^M^Mf^ = 8(l- e ) 

spins 

The one-loop exchange contribution, evaluated to all orders in e, is given by 

£ M<pM<f>1 = £ MfM^+A^ , (2.33) 



u 2 



(2.32) 



spins 



spins 
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where 
id) 



1 - 2e)—( (1 - 4e + e 2 )t 2 - 2e(2 - e)tu + (1 - e) V^Box^ (s,t) 



+ 8(1 - 2e)^(e(2 - 3e - e 2 )t 2 + 2e(l - 3e - e 2 )tu - (2 - 2e + 3e 2 + e 3 )u 2 )Box (6) (s, u) 
2e(l-e)(u 2 + esi)iV / 



s 



(1 - 2e)(3 - 2e) i 



(3 - 2e) (2e(l + e 2 )t 2 + e(3 + 2e 2 )tu - 2(1 - e + e 2 )u 2 ) 



Tri(s) 



8 1 

+ ( e(2 - 5e + 2e 2 - e 3 )i 2 + e(l - 3e + e 2 - e 3 )tu - (1 - e)(2 - 3e - e 2 )u 2N )Tri(t) 

1 — 2e s V / 

- T A r -(e(2 - 4e + e 2 - e 3 )t 2 + e(2 - 3e - e 3 )tu - (1 - e)(2 - 4e - e 2 )u 2 )Tri(u) . 

(2.34) 

Using these results, and the computation of the two-loop exchange terms, we again find that the 



additional singular terms in Bhabha scattering are described by eq. (|2.9| ), where (not surprisingly) 



Hrs is given by precisely the same expression ( 2.21 ) that we found for e + e 



MA* 



2.5 Finite Contributions to the Amplitudes 
2.5.1 e + e~ -» fi+jT 

Finally we give the real (dispersive) part of the finite remainder in eq. (^ 
amplitude in the CDR scheme. First we treat the e + e~ — ► process (2.14). It is convenient 

to decompose the finite part according to the number of light flavors, Nf, 



interfered with the tree 



J2 Re 



M 2)fin M 0)t 



F (0) + Nf F (l) + N 2 F (2) 



(2.35) 



spins 



In the s-channel, the functions are given by 



f® + s\f®' 



(2.36) 
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where 



F (o) = 2 ^X 2 + (x 2 + y 

y 



4(2 - 0(Lh(-x) - XU 3 (-x) + ^X 2 U 2 (-x)) + ^vr 2 Li 2 (-x) 
+ (hi + OX 3 + \{2- OX 2 Y - l -XY 2 - 3(X - Y)X - 2^((3 + $X - 3Y 



+ I(n - 16£)X - %Y + ef-12C 3 + tt 2 + ^ 



__ 43 A 15. 29 2 511 



+ (x -y) 



8Li 4 (-x/y) + 6(2 + QLU(-x) + (4(1 - £)X - 12F)Li 3 (-x) 



+ 



(2 - OX 2 - 4XY + -tt 2 jU 2 (-x) - (2 + (Lia(-x) - XLi 2 (-: 
-1(6 + OX 3 + \x 2 Y + 1(1 + 4^)vr 2 X + 1(10 + t)X 2 - 1(2 - £)XY 



1 



7T 



+ -(1 + 6£)* - 4(2 - £)Cs - (1 + 40— - 6£ U + £ -6C4 - 2Cs + 2 



7T 



+ (2 



Li 3 (-x)-XLi 2 (-x)) + (-(5-3e)^--(3-e)Xy-(l-40y )X 



+ 



1 



(1 + 60*- iy-6£)X-4C 3 + 



7T - 



3 ' 



(2.37) 



9' 



{x 2 + y 2 " 



) f £ [l2(Li 3 (-x) - XLi 2 (-x)) + (4X 2 - 6X(Y + 2 lnfoi 3 /*)) + 20^ 2 



87 



35 



29X + 361n(/i 2 /s) +33 )X + — \n(^/s) + — Cs + -riT + 



685 \ 



2£(x - y) {^(X 2 - 3(X - 1) ln(/i 2 /s) - yl + 4^ 2 + 8)x - 2tt< 
£(3X + 61n(/i 2 /s) + 16)xl , 



F (2) 

with 



lnV/-) + ylB(M 2 /-)-^ + f 



x = - , y=~, X = ml -- 



s s 

and the symmetry operation S is given in eq. ( |2,25 ) 
In the u-channel, the are given by 



r = i„|-||, 



(2.38) 
(2.39) 

(2.40) 



_p(i) _ p(i) 



(2.41) 
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2 X —y ({V - Wf + vr 2 ) - 2(x - y) Q - 3?) V 

+ (a; 2 + y 2 



(2 + (Lu(-v) - VLi 3 (-v) + ^ 2 Li 2 (- 
+ (2 - ^Li^-v/u;) + (V — W)(Li 3 (-v) + Li 3 (-™) - WU 2 (-w) 

-hv + w)u 2 (-v)) 



+ 2(6 + 0yLi 2 (-t;) 



1 



(x-y) 



- \^V Z W + (4 - £)^ 2 ^ 2 - \{7 - 2i)VW' 4 - ^(1 - 2£)^ 4 

1 93 

+ 6Vi^ 2 - 3VF 3 + i6^yvF + -(9 - 160 w 2 + -^i w 

+ vr 2 (-2V 2 + |(3 + j)yw - 1(3 - £)1^ 2 - 6V - (3 - £)W 

/ \ 7T 2 ^11 

+ 4((2 - 0^ - 2(1 + OW) Cs + 34C 4 - 15Cs - (25 + 96£)^ + 

6((2 + £)Lu(-v/w) - (2 - £)Li 4 (-w)) - 16Li 4 (-w) 

+ 4((1 - OW + (2 + 0^)Lia(-«;) + i^iV - (2 + £)^)Li 3 (-?;) 

- (2 + e)(Li 3 (-u») + WT,i 2 (-w)) - 4(Lia(-«) - VLi 2 (-t;)) 

+ (-4V 2 + 2(2 + £)VW + (2 - £)1^ 2 - (10 + 3^)y^Li 2 (-«) 

+ ^V 4 - |(2 + 0^ 3 ^ + ^(2 + £)^ 2 - ^(4 + 70VW 3 + i(2 + £)1^ 4 

+ -£^ 3 - (2 + ^)y 2 VF + 2VW 2 - -(10 + £)VF 3 + (5 + 6£)^W - -(5 + 6£)W 2 
3 6 2 

2 

+ y (2(6 - i)V 2 - 2(13 + 4£)VW + (5 + 3£)W 2 + 2(12 - £)V - (3 - 2£)w) 

2 

- 4((2 - 0^ - w)Cs - 6^(2^ - + (121 + 3£)C4 + ( 2 + 5 0Ca - (15 + 26£)^ 
+ 2^(Li 3 (-i>) - VU 2 (-v) - V 2 W) - (2 - 0(Li3(-w) - WU 2 (-w)) 

- (1 - 0VW 2 + ^(5 - + 6£,VW + hi - 6£)W 2 - 6£W 

6 2 

2 2 

+ y (2(4 + 30^ + (8 - 5£)W) - (6 + OCa + (1 - 18£)y , 

(2.42) 
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5 (* 2 + i/ 2 ) 



121 Li 3 (-u;) - TyLi 2 (-w) + 2(Li 3 (-v) - FLi 2 (-t;)) 

- -V 2 W + -Fiy 2 - -W 3 + {W 2 - 2VW -3W + vr 2 ) ln(/i 2 /(-s)) 
2 2 3 

+ 29(W - 2V)W - 33iy + 2tt 2 (3V - AW) 



tt 2 1370 



+ 871h(m7(-«)) + (35 + 12£)C3 + (7 - 58£) — + () 

2 y 3 + 2(V - Wf + 13(2(1/ - W)V + IF 2 ) 
+ 6(2(V - + VF 2 + 2F - + vr 2 ) ln(/i 2 /(-s)) + 16(2F - W) 
-2tt 2 (V + W) + 9tt 2 

(3W -6V + 6 ln(fi 2 /(-s)) + 16) W + 3tt 2 



Fi 2) = ^(x 2 + y 2 ) 



lnV/(- S )) + yln(^/HO) + f 



(2.43) 
(2.44) 



Here x, y are defined in eq. (|2.4C| ), whereas v, w, V, W are denned in eq. ( 2.27 ). 

In the t-channel, the functions are given by the action of the symmetry S of eq. ( 2.25| ) on 
the u-channel results, 



fW = s !>« 



5=1 



(2.45) 



The two-loop virtual contribution to the e + e — > unpolarized cross section, restoring 

overall factors and averaging over initial spins, is given by 



da® 
dt 



1 (Aira) 
x 



167TS 2 



'M a) M 0)t 



(2.46) 



spins 



2.5.2 Bhabha Scattering 

For the finite two- loop remainder for the Bhabha scattering process ( 2.29Q , we quote only the 



(s <-> t) symmetric sum of the two exchange terms required by eq. ( 2.30 ). Again we decompose the 
answer according to Nf, 



£ Re 



spins 



+ U 



Re 



F (0) + Nf F W + N 2 p{2) 



(2.47) 



In the s-channel, the functions FW are given by 

p(i) _ 



(2.48) 
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where 



2^Y 2 - 2x z ({X - Yf + vr 1 



x- 



+ 



X 



-4(Li 4 (-x/y) - LU(-y) + XLi 3 (-y)) + 2(4Y - 3X - l)Li 3 (-x) 

+ A(X 2 - 2XY + \x + vr 2 )Li 2 (-x) + \x A + \x 3 Y - 4X 2 Y 2 + \xY 3 - \y a 

23 3 93 

X 3 + -X 2 Y + 9XY 2 - 6F 3 - 5X 2 - 21XY + 23Y 2 + — (X - 2Y) 

12 2 i y ' 

vr 2 

+ y (-17X 2 + 32XY - 18Y 2 - 17X - 26Y) 



2C 3 (3X - 8Y) + 15C 4 - 38C 3 + 47^- + ^ 

o 8 



+ 



- x) 



12 



-lOLU(-x) + 6XLi 3 (-x) - ( X 1 + 2yJLi 2 (-x) + — X 



+ y " 2 ((x " y) + } + 2 + r 6C3 + 2 + 12 ) x + 20a 



+ 16(Li 4 (-x/y) - Li 4 (-y) + XLi 3 (-y)) + 4(3X - 2Y - 2)Li 3 (-x 

- 4(X 2 - 2XY -2X + vr 2 )Li 2 (-x 
5 



2 o o \r i _2\y ; / ™N — X 3 Y + 8X 2 Y 2 — XY 3 + —Y^ 

12 3 3 3 



+ -X 3 + X 2 Y + 18X 2 - 4XY + 2Y 2 + — (11X 2 - 20XY + 4Y 2 + 9X) 
6 3 

4( 3 (3X - 2Y) + 88C 4 + 8Cs + 2vr 2 , 



(2.49) 



9 1 x 



36(Li 3 (-x) -XLi 2 (-x) - -(X-4)Xln(/i 2 /s) 



+ i(y - 3)y(ln( /U 2 / S ) + ln( M 2 /(-t)))) + 2X 3 - 24(x + ln(/i 2 /(-t)))XY 



+ 12XY^ - 8Y 3 - 19X^ - 58(X - Y)Y + tT 18X - 28Y + 12 ln(^/(-t 



+ X - 66Y + 87(ln(/i 2 /s) + hu> 2 /(-i))) + 58C 3 + 44vr 2 + 



2740 



- 2x 
-4y 



3(ln(/i 2 /s) + ln(^ 2 /(-t))) + 16JX 2 - 6vr 2 X 
X 3 + hln{fi 2 /(-t)) + 8)x 2 + f3(ln(/i 2 /s) + ln(/i 2 /(-t))) - 2tt 2 + 16 )x - 3ir 2 



p{2) = W 
9 x 



(2.50) 
(2.51) 



lnV/ S ) + InV/H)) + f (intf/s) + ln(M 2 /H))) " - 2 + y 
and x, y, X, Y are denned in eq. d2.40[) . 

In the i-channel, the functions are given by the action of the symmetry U of eq. ( 2.31| ) on 
the s-channel results, 



(2.52) 
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In the u-channel, the functions are given by 

p(€) = pW +u [p® 

where 



(2.53) 



pW 



F^ 



y 2 
+ — 

X 



-4(m(-t;) + WLi 3 {-v)) + 6V(Li 3 (-v) + Li 3 (-w)) + 2Li 3 (- 
+ (~~V 2 - WW + ^W 2 -V + W + lly^Li 2 (-v) 



13, 



11 23 21 

+ ( i y3 ~ n y2w ~ ^r vw2 ~ 7^ y2 + ~r vw -5V + —W - 12Cs + — 



12 



31. 

T 



93 



93 

T 



47 



vr 2 511 



V + — C4 - 20C 3 - 109— + 

8 12 16 



+ 



?/(l - x) 



x 



5Lu(-v/w) + 6V (Li 3 (-u) + LiaC-u;)) 

- I ($V + W- 2){V -W) + b^-]u 2 {-v 



+ (_V + -V 2 W + -V 2 - -VW + -V-12 + — (V + 7)) V 
V 4 2 12 4 2 12 v 7 

+ 16(Li4(-«) + WLig(-w)) - 12V(Lig(-t;) + Li 3 (— «;)) + 8Li 3 (-w) 
+ f(3V + W - 4)(V -W) + h^]u 2 (-v) + ( — —V 3 + 2V 2 W - %VW 2 + - V 2 + -VW 



TT 



+ 18V - 161V + — (-4V + 9TV - 7) V - — C4 + , 



12 
253 



1 |y_ 

9\ x 

— x 6 



2y 



2 r 



6 



36(Li3(-«) - VLi 2 (-u) ) + 871n(/x 2 /(-s)) + 47Ca + 



(2.54) 



1370 



((V - W) 2 + tt 2 + 2{V - W)) ln(/x 2 /(-s)) 

+ (13V - 31JV + IOtt 2 + 65) (V - W) + 16vr 2 

-2V 3 + 18V 2 1V + 12VTVln(/i 2 /(-s)) + 9V 2 + 20VW + 18(V + W) ln(// 2 /(-s)) 
2 



+ 33V + y( 4V ~ 121n (^ 2 /(-«)) - 29 



^[ln 2 ( M 2 /(- s) ) + ^ln(^ 2 /(- S )) + | 



(2.55) 
(2.56) 



Here x, y are defined in eq. (|2.4C| ), whereas v, w, V, W are defined in eq. ( 2.27Q . 
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2.6 Checks on the Result 



We performed several checks on our calculation. The calculation was performed with the computer 
algebra programs Maple, Mathematica, and FORM. To check the code, large parts of the calculation 
were performed independently with alternative programs written in different languages. Various 
checks were applied to the integral reduction procedures described in refs. |l4|, [l^, 18, |l9| and 

our implementation of them. For example, we reproduced the double box ultraviolet divergences 
in d = 8 and d = 10 reported in ref. ||, and several other previously calculated double box 
integrals fl(| . An additional check on the non-planar tensor integrals is that unphysical l/(t — u) 
poles occur in the representation of these integrals in terms of the master integral basis we used ]l5| ; 
however, in the series expansion in e such poles drop out after delicate cancellations between the 
various terms. 

We checked the gauge invariance of the scattering amplitude by explicitly calculating the Feyn- 
man diagrams in a general £ gauge and observing that the gauge dependence drops out in the 
final result. This provides a non-trivial check of the diagrams and parts of the integral reduction 
procedure. 

A strong check on the final result is provided by the matching of the IR divergence structure 



of the two-loop scattering amplitude with Catani's formula (|2.9|) , as discussed in sect. |2.l| . A 
given integral will contribute to both infrared divergences and to finite terms. Thus a check of the 
divergent terms provides an indirect check that the finite terms have been correctly assembled. 

Finally, we observed for small scattering angles a suppression of the leading logarithms, I = 
ln(# 2 /4), e.g. in the limit s — > in the i-channel for process ( |2.14j ). In other small-angle limits 
(those not enhanced by the photon propagator pole) the leading power-law behavior is of course 
less singular, but it is dressed by large logarithms of the type and Nf£ 3 . But in the t-channel 
s — > limit it cancels down to £ 2 and Nrl. This behavior is in accord with a generalized eikonal 
representation for small-angle scattering [0]. 



3 Conclusions 

In this paper we presented the two-loop QED corrections to e + e~ — > and to Bhabha scat- 

tering. We presented the results in terms of two-loop amplitudes interfered with tree amplitudes 
and summed over spins in the context of conventional dimensional regularization. In these results 
we have set the small electron and muon masses to vanish. (This is an excellent approximation for 
the highest energy current and future electron-positron colliders.) 

The two-loop amplitudes presented in this paper are infrared divergent. To make use of them in a 
Monte Carlo program for the NNLO terms in the cross section, they must be combined with lower- 
loop matrix elements including photon emission, which should be computed using conventional 
dimensional regularization, at least in the singular regions of phase space. In particular, the pieces 
that need to be computed (for the Bhabha case) are 
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• the e + e~e + e~ one-loop amplitude interfered with itself. 

• the e + e _ e + e~7 one-loop amplitude interfered with a five-point tree amplitude, and 

• the e + e _ e + e~77 tree-level squared matrix element, 

The interference of the dimensionally regularized one-loop four-point amplitude with itself does 
not appear to be in the literature. Nevertheless, it should be relatively straightforward to obtain, 
given that it involves only one-loop amplitudes with four-point kinematics. The required integrals 
are given to sufficiently high order in e in eq. ( [2.26 ). 



The QED one-loop five-point amplitude interfered with the five-point tree is a rather involved 
object to compute from scratch. However, the closely related one-loop helicity amplitudes for one 
photon and two quark pairs are known []33| , [34| , |35f | , and it is a relatively simple matter to modify the 
color factors to obtain the corresponding QED amplitudes. The one-loop helicity amplitudes are 
in the 't Hooft-Veltman scheme. They can be converted to conventional dimensional regularization 
by altering the tree amplitude appearing in the coefficient of their singular terms |36|]. Thus the 
e + e~/x + /i~7 and e + e~e + e~7 one-loop amplitudes may be extracted from the known literature 
through O(e ). 

Because of the 1/e 2 infrared divergences that are encountered in the phase-space integral, in 
regions where the photon is soft or collinear, one might seem to require the one-loop five-point 



amplitude through 0(e 2 ). However, this is not necessary [21|. Instead, one can replace the five- 
point amplitudes in singular phase-space regions by a combination of four-point amplitudes (which 
are given in this paper to the required order in the dimensional regularization parameter) and 



splitting amplitudes |57], p8 |. The one-loop splitting amplitudes for QCD are enumerated to the 
required order in refs. |21]; the case of QED follows as usual by an appropriate conversion of color 
factors. 

The tree-level helicity amplitudes for e + e~ and e + e _ e + e~77 have been known for a 



while p9fl . (They also can be converted from the four-quark two photon amplitudes in ref. [35 
for example.) In infrared-divergent regions of phase space one must include higher order in e 
contributions from the matrix elements. Systematic discussion of these regions, where two particles 
can be soft or three collinear, has been presented in refs. [^] for the case of QCD. Once again the 
results for QED can be obtained by a conversion of the color factors. 

Even with all of these matrix element ingredients assembled, it is a nontrivial task to devise a 
numerically stable method for carrying out the singular phase-space integrations. Nevertheless, this 
task is very analogous to that required to obtain QCD jet predictions at next-to-next-to-leading 
order, so it is likely that it will be attacked soon. 

Besides the obvious application of the present paper to refined theoretical predictions for Bhabha 
scattering and for electron-positron annihilation into muons, it also serves as a further test of 
methods that can be applied to analogous QCD processes. We are confident that many more 
multi-particle two-loop amplitudes will be calculated before long. 
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